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Abstract 

Using liquid 3 He in aerogel as an example, it is shown that correla- 
tions in positions of impurities affect the temperature T c of transition 
of Fermi liquid in an unconventional superfluid or superconductive 
state. The effect is significant if the correlation length for impurities 
is greater than the coherence length in the superfluid or superconduc- 
tive state £o- F° r 3 He in aerogel the suppression of T c is expressed 
in terms of the structure factor of aerogel. With the account of the 
fractal structure of aerogel a simple expression is obtained for the de- 
crease of T c from its clean value. This expression is in a satisfactory 
agreement with the experimental data. 

1. For unconventional Cooper pairing the temperature of transition into 
the superfluid (or superconductive) state is suppressed both by magnetic 
and by non- magnetic impurities pQ. The amount of suppression in most 
cases is well described by the theory of superconductive alloys of Abrikosov 
and Gor'kov (AG in what follows) [2] . Superfluid 3 He is the best understood 
example of unconventional Cooper pairing. The Cooper pairs here have 
orbital momentum 1 = 1, and spin s = 1. Floating impurities can not 
be introduced in the liquid 3 He - they stick to the walls. To get around 
the difficulty a high-porosity silica aerogel has been used as an impurity 
[3j H] . The aerogel is a self-supporting network of silica strands with a typical 
thickness 3-4 nm. Rigidity of the network assumes existence of correlations 
in positions of its elements. The homogeneous scattering model (HSM) [3], 
which is a generalization of the AG theory for the p-wave pairing, does not 
take into account these correlations. This drawback is essential, since HSM 
does not provide a quantitative description of properties of superfluid phases 
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of 3 He in a presence of aerogel. In particular, HSM does not describe correctly 
the dependence on pressure of the magnitude of suppression of superfluid 
transition temperature of 3 He in aerogel T a with respect to its bulk transition 
temperature T b . 

In the AG-theory, as well as in the HSM, the relative value of the sup- 
pression r ba = Tb ~ Ta is determined by the only parameter x = £o/kr, where 
£o = hvF/(2TrT a ) is the superfluid coherence length and l tr is a transport 
mean free path. For r ba <C 1 the theory predicts r ba = TiT - Both £ and 
l tr can be found from independent experiments [6]. The measured transition 
temperature T a is greater than the value, calculated within the HSM with 
the use of the known values of £ an d h r - The other suggested models [5J [7j 
invoke an effect of restricted geometry as a mechanism of suppression of T a . 
One of the models considers 3 He in a gap between two diffusely scattering 
planes, and the other - 3 He in a spherical void. Although a better agreement 
with experiment can be achieved in this way, relation of the models to the 
real aerogel remains unclear. 

In a present paper effect of correlation in position of elements, forming 
aerogel, on the superfluid transition temperature of 3 He is considered directly. 
The effect is the stronger the greater is the ratio of correlation length in 
aerogel R to £ - When R 2 ~ l tr £ the change of r ba stemming from the 
correlations can be of the order of the original r ba . This region is of particular 
interest for the aerogels, used in experiments. To simplify the argument we 
consider here only a region r ba -C 1. 

2. At the p-wave Cooper pairing the order parameter is 3x3 complex ma- 
trix A^j. The first index assumes tree values enumerating three projections 
of spin of the Cooper pair. The second index enumerates three projections 
of its orbital momentum. To find T a assuming r ba <C 1 one can use Ginzburg 
and Landau equation: 

-r A()+W ,)-|e(^ + 2^)=0. (1) 

Here r = (T b — T)/T b <C 1. Effect of impurities (aerogel) is introduced via a 
real symmetric tensor %(r), it describes a local depression of the transition 
temperature and its possible splitting for different orbital components of 
Apj. Interaction with the impurities can be written in Eq.(l) in a local form 
because tensor %(r) varies on a distance ~ £o ; while A^ in a vicinity of 
T b varies on a distance ~ £(T) = £o/v^" ^ £o- A coefficient in front of 
the derivatives is written as in the Ref. [8], i.e. £f = -^7^£o — 0,7^- At 
t > Eq.(l) can have physically meaningful solutions. These solutions are 
generally speaking nonuniform and can be localized. By the definition T a is a 
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temperature of the onset of the long-range order in liquid 3 He. In the present 
formulation it corresponds to the minimum r = r a , at which Eq.(l) can have 
delocalized solution. The equations (1) for different spin components /i are 
not coupled and the spin index is not relevant. A problem of finding r a b turns 
out to be analogous to a problem of finding mobility edge of a spin-1 particle 
in a random potential % (r), and r is analog of energy. To avoid cumbersome 
calculations let us introduce one more simplification. Instead of the sum of 
two gradient terms in Eq.(l) we use a model isotropic expression with only 
one term and an "average" coefficient 



rA llj + A lim {T)-g[-=rf)=0. (2) 




The average £f is defined in the following way. For the longitudinal compo- 



nent of A^i, which meets the condition = 0, the coefficient in front of 
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the derivative in Eq.(l) is (3/5)£f. For two transverse components, defined 
as A^ = where ip^ are scalars, a proper coefficient is (9/5)^. The aver- 
aging over three possibilities renders £f = (7/5)£f. Since ^ ~ 0, 7£q with a 
reasonable accuracy £| ~ £q- In what follows £q will be used as a coefficient 
in front of the gradient term in Eq. (2) . 

Explicit form of tensor rjij (r) depends on the particular structure of aero- 
gel. We assume here, that aerogel consists of spheres of uniform radii p, 
distributed with an average density n. The values of p and n are chosen to 
provide the required values of porosity and of the mean free path for single- 
particle excitations. It will became clear later, that a concrete form of the 
elements is not essential, but for explicit calculations spherical form is prefer- 
able. For the porosity 98,3% and k r ~130 nm p «2 nm, i.e. p <C £o- For 
such model 

= Y, f lj 1 i\r-r s ), (3) 

s 

where r s is a coordinate of the sphere number s, and T)j]\r) a potential 
induced by a single sphere, it depends on the cross-section of scattering of 
the single-particle excitations on the sphere [9]. Assuming for the sake of 
definiteness diffuse scattering, we arrive at r ^> p at the expression [10J: 



(i) / \ P ~ ~ i 



tanh 



(4) 



According to Eq.(4) rff (r) decays on a distance ~ £ - At r < £ Vjii'-j ~ 
(p/£o) 2 - Let us treat ^(r) as a perturbation. Green function for Eq.(2) after 
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Figure 1: 



averaging over realizations of 77^ (r) has a form: (G mn (T; k, k')) = (27r) 3 5(k — 
k')5 mn G(T; k), where 

The self-energy S(r, k) up to the factor (27r) 3 <5(k — k')5 mn is the averaged 
sum of the series shown diagrammatically in Fig.l. The arrows in the figure 
correspond to the unperturbed Green functions 

^(r;k) = -^. (6) 

As usual, the integration over momenta of internal lines is assumed. Wavy 
lines correspond to the Fourier transform of the potential rjji(r): 

^(k-kO^^^k-kOE^'-^ 3 , (7) 

s 

where k and k' - momenta, corresponding to the arrows coming in and out 
of the vertex. The averaging is performed over coordinates r s of particles 
forming aerogel. The "mobility edge" is found as a pole of the Green function 
Eq.(5) at k = 0: 

t = E(t,0). (8) 

Fourier transform of r)j\\~k — k') can be found directly from Eq.(4). Here we 
need only r]fi (k — > 0) = 5jii]^(0), where: 

^ (1) (0) = yP 2 6. (9) 
Consider consecutive terms of the series Fig.l. The first order term is: 

( %i (k-k')}=#(k-k')£e^). (10) 



Assuming that aerogel is on the average uniform we have: 



(]T e*^-^'*) = (27r) 3 (5(k - k>. (11) 

s 

As a result, in the first order on the perturbation 

i 1) =W 1) (0) = ^- (12) 

For the comparison with HSM the answer is expressed in terms of the trans- 
port mean free path of the excitations h r , which is defined as: 

r = ~Vn. (13) 

This definition corresponds to a diffuse scattering of quasi-particles by the 
randomly distributed uniform spheres with radii p. The first order correc- 
tion (12) coincides with the result of HSM for small Tf, . The second-order 
correction is: 



d 3 h 

V s ' 

(14) 

Instead of the unperturbed Green function Eq.(6) the average Green function 
Eq.(5) is substituted in the r.h.s. of Eq.(14). This substitution is known as 
the self-consistent Born approximation [TT]. The summation in Eq. (14) is 
over the relative coordinate r st = r s — r t . The averaged sum is the structure 
factor: 

(£e Kkl ~ k)rst )=S(ki-k), (15) 
t 

it characterizes correlations in positions of the particles, forming aerogel. A 
structure factor is directly measured by the small-angle x-ray scattering. For 
the isotropic aerogel S(k) does not depend on a direction of k. Substituting 
k = in Eq. (14) we arrive at: 

T ba^ji = n J ^m(-ki)^2(ki)5(/c 1 )G(r 6a ,A;i)^^. (16) 

Following Eq. (8) and neglecting a possible small change of the spectrum we 
arrive at G(Tb a , ki) = — l/(£ ^i) 2 - At evaluation of the integral in Eq.(16) an 
account must be taken of the fact that in an interval of scales g < r < R, 
where R is an upper limit of the interval, aerogel has fractal structure. It 
means, that in the corresponding interval of wave vectors k m i n < k < k max its 
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structure factor S has a power dependence on k: S ~ k D . The exponent D 
is referred as a fractal dimension. In particular for the samples of 98% aerogel 
used in Ref.p2] k min -5^7- lCT 3 !" 1 , k max ~ 1 2 • lO^A -1 , D«l,8. 
This rate of decay of S'(fc) secures convergence of the integral in the infinity. 
A principal contribution to the integral is provided by the region of small k. 
That allows to substitute in the integrand rf^J}a — ^i) — Vjmfi) = SjiV^ifyi 
then 

i 2) = ^ (1) (0)] 2 I S(h)G(0, (17) 



(2tt) 3 ' 

Therefore a particular shape of the elements, forming aerogel is not essential. 
Within the fractal interval the integral in Eq.(17) is ~ / dk/k D and at k min — > 
it diverges at small k. Actually the integral converges since at k <c k m ^ n 
dependence S(k) saturates. That reflects an absence of correlations on a 
distances exceeding R ~ l/k min . A smooth cut-off can be introduced by 
assuming a simple model law of decay of correlations, depending on R. Then, 
by definition, R is a correlation radius. The structure factor 5*(k) can be 
expressed in terms of the pair correlation function C(r): 



S(k) = (5> lkrst ) = n [ C(r)e ikr d 3 r. 



At r <C R for a fractal with the dimension D correlation function behaves as: 
C(r) fa A(R/r) 3 ~ D , where A - is a coefficient. At r ^> R correlations vanish 
and C(r) 1, i.e. the integral in Eq.(18) diverges. It converges if C(r) is 
substituted by v(r) — C(r) — 1. This substitution corresponds to separation 
of the effect of correlations from that of uncorrelated distribution of elements. 
The unity gives a contribution to S'(k), which is proportional to S(k) and 
does not affect the following calculations. The integral in Eq. (17), denoted 
as Q in what follows, with the aid of straightforward transformations can be 
expressed in terms of v (r): 

Q = J S{k l )G{k l )^ = -»/ v(r)rdr. (19) 

At r R function v(r) ~ A(R/r) 3 ~ D . For D > 1 the integral in Eq.(19) 
converges at small r, therefore the fractal asymptotic can be used up to r = 0. 
At r > R vanishing of correlations has to be taken into account. Following 
Ref. P2] we assume here the exponential decay of correlations, i.e. substitute 
for v(r): 

v (r) = [A - - l]exp(-r/R). (20) 
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Coefficient A is found from the normalization condition: 4nn J v(r)r 2 dr = 
-1. For aerogels in question nR 3 > 1. Then A = 2/T(D), where T(D) 
is Euler Gamma-function. With this A we arrive at: / v(r)rdr = R 2 (3 — 
D)/{D - 1). As a result 

Suppression of the transition temperature with the account of the second 
order term 

Tba ~ 4kr ^kr{^0krD-l) ^ 

turns out to be smaller then that, predicted by HSM. The factor 

Ti^rrfer S rows when £ decreases. Compensation of the small parameter 
Co/hr by a big ratio (i?/£o) 2 can result in a product of the order of unity. In 
this case higher order terms in the expansion of Tb a over rjji have to be taken 
into account as well. The third order term is: 

i 3) = ^ (1) (0)) 3 / ^ ^ 2 r G(h)G(k 2 )(j: e *^ e *^). (23) 

The averaged sum here depends on three-particle correlations. The higher 
order terms respectively depend on the higher order correlation functions. 
The sum of the series can be found if the assumption is made, that the higher 
order correlation function can be decoupled in products of the two-particle 
correlation functions. In particular: 

e ^ r -e ik2r " s ) = 5(k 1 )5(k 2 ). (24) 

t,u 

In this case rf^ is a product of r^f by r]^(0)Q etc.. The consecutive terms 
form geometric series with the sum: 

r (i) 

Tba = T^Twq- (25) 

With the given above expressions for and Q we arrive at 

^4 Zohr D-l 

The obtained expression for Tj, a in both limits of large and small R gives 
physically natural results. At R <C £o the result of HSM is reproduced. If, 
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on the other hand, R is so large that f 3 - 

Hr 




3> 1, then the mean free path 



l tr drops out of the expression for the transition temperature: u a ~ 



Such behavior is in line with the qualitative argument, based on the presence 
in aerogel of the low-density regions, or 'voids', with a characteristic size £ a 
[3 E] . According to the argument for £ < £ a superfluidity sets on starting 
from interior of the 'voids', then r ba ~ (£o/£a) 2 , while at £ 3> £ a homogeneous 
limit is recovered T ba ~ ^o/hr- For interpolation between the two limits Sauls 
and Sharma [H] suggested to substitute in the formula of HSM for T a instead 
of the pairbreaking parameter x = ^o/2l tr an effective pairbreaking parameter 
x — x/(l + Ca/ X ) with ( a = ^ a /ltr- The obtained expression turns out to be 
in a good agreement with the data for 98% aerogel. Eq.(26) can also be 
rewritten as the HSM formula (12) with a substitution of the parameter 



x instead of x. To do this one has to set £ a = Rny(3 — D)/8(D — 1) in 
the definition of x. Therefore Eq.(26) with a proper choice of R has also 
agree with the data for T{, a , when is small. The constraint is not crucial. 
The method of correlation function [TS] makes possible generalization of the 
obtained result for finite t^. 

Good agreement with the data for T a in a 98% aerogel p3] is achieved 
at £ a ^500A and k r ~1400A. With these values of parameters at pressure 
above ~20 bar C^/x ~1,6, i.e. effect of correlation is essential. 

In conclusion, one can see that the account of the correlations in positions 
of the elements forming aerogel explains suppression of the temperature of 
transition of 3 He in the superfluid state by aerogel. The observed difference 
in properties of superfluid 3 He in aerogel from the predictions, based on the 
HSM [16], indicates that correlations are essential for a proper interpretation 
of these properties as well. 

The above argument can be applied to unconventional superconductors 
as well. In particular, for some of high-T c superconductors the AG-theory 
does not properly describe the suppression of T c by impurities. The nu- 
merical analysis of the Ref . [17] demonstrates that the discrepancy originates 
from the fact that the correlation length £ in these materials is significantly 
smaller than in conventional superconductors and can be comparable with 
the average distance between the impurities. 
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